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ABSTRACT 

We describe observations of the 2003 outburst of the hard-spectrum X-ray transient 4U 1901+03 
with the Rossi X-ray Timing Explorer. The outburst was first detected in 2003 February by the All- 
Sky Monitor, and reached a peak 2.5-25 keV flux of 8 x 10 -9 ergs cm -2 s -1 (around 240 mCrab). The 
only other known outburst occurred 32.2 yr earlier, likely the longest presently known recurrence time 
for any X-ray transient. Proportional Counter Array (PC A) observations over the 5-month duration 
of the 2003 outburst revealed a 2.763 s pulsar in a 22.58 d orbit. The detection of pulsations down 



ergs cm s (2.5-25 keV), along with the inferred long-term accretion rate of 



to a flux of 3 x 10" 

8.1 x 10" 11 Mq yr -1 (assuming a distance of 10 kpc) suggests that the surface magnetic field strength 
is below ~ 5 x 10 11 G. The corresponding cyclotron energy is thus below 4 keV, consistent with the 
non-detection of resonance features at high energies. Although we could not unambiguously identify 
the optical counterpart, the lack of a bright IR candidate within the 1' RXTE error circle rules out 
a supergiant mass donor. The neutron star in 4U 1901+03 probably accretes from the wind of a 
main-sequence O-B star, like most other high-mass binary X-ray pulsars. The almost circular orbit 
(e = 0.036) confirms the system's membership in a growing class of wide, low-eccentricity systems in 
which the neutron stars may have received much smaller kicks as a result of their natal supernova 
explosions. 

Subject headings: accretion — pulsars: general — pulsars: individual (4U 1901+03) — X-rays: stars 



1. INTRODUCTION 

More than half of the known X-ray pulsars are tran- 
sient sources which were discovered during bright out- 
bursts. Most of these are Be-star/X-ray binaries, which 
are believed to be progenitors of double neutron st ar bi- 
naries (e.g. IBhattacharva fe van den Heuvell|l99H) . In 
contrast to binaries with low-mass (< 1 Mq) compan- 
ions, the neutron stars in these high-mass X-ray bina- 
ries (HMXBs) typically accrete through the companion's 
stellar wind. The mass transfer rate M — as well as its 
variability — depends in a sensitive manner on the wind 
properties, as well as on the neutron star spin period 
and magnetic field strength. A smaller class of pulsars 
accrete from supergiant companions, which may fill their 
Roche lobes and thus accrete persistently. Tidal forces 
in these binaries act to circularize the orbits on a time 
scale much shorter than the active lifetime, while the 
wider, wind-accreting bin aries typically have m oderate 
to high eccentricities (e.s. iBildsten et aTlll997j) . In re- 
cent years, a third class of binaries has emerged, with 
wide > 20 d orbits but low eccentricities e < 0.1. These 
sources present difficulties for the commonly accepted 
formation scenario in which the natal supernova event 
imparts a "kick" to the neutron star, leading to an ini- 
tially eccentric orbit. The tidal forces which act so ef- 
ficiently in the Roche-lobe filling systems cannot circu- 
larize wider orbits within the source lifetime, suggesting 
that the initial kick s in these wide, cir cular binaries may 
be unusually small ijPfahl et alJl2002f) . 

For some transients the interval between outbursts can 
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be as long as 20 yr. Little is generally known about 
sources with such long duty cycles, due to the dearth 
of observations (in particular with large-area modern 
instruments with good timing capabilities). For some 
candidate X-ray pulsars no pulsations have even been 
detected, and the classification comes from a hard X- 
ray spectrum, typical in confirmed pulsars. 4U 1901+03 
(I = 37?16, b = — 1?25) is such a source, previously de- 
tected just o nce before in outburs t by Uhuru and Vela 5B 
in 19 70-1 l|Forman et alJ Il976t iPriedhorskv fe Terrelll 
Il984[) . Due to the hard spectrum measured during those 
observations, the source was tentatively identified as an 
HXMB. Consequently, we selected this source as one of 
a group of hard transients with positions known to 10' 
or better, as candidates for target-of-opportunity obser- 
vations by the Rossi X-ray Timing Explorer (RXTE). 
In 2003 February a new outburst of 4U 1901+03 was 
detected by the All- Sky Monitor (ASM) aboard RXTE 
l)Oallowav et alJl2003b(i . The source was also detected 
with the IBIS and JEM-X hard X-ray instruments 
aboard I NTEGRAL between 2003 March 10 and April 
13 2003 l|Molkov et all 12001 . An RXTE Proportional 
Counter Array (PCA) scan across the Uhuru position 
led to more precise coordinates of R.A. — 19 03 m 37!l, 
deck = +3°11'31" (J2000 . 0), with an estimated 90% 
confidence uncertainty of 1' l|Gallowav et all2f )03a) . Fol- 
lowup pointed RXTE observations detected coherent 
pulsations with a period of 2.763 s. Variations in the 
observed pulse frequency were also observed, suggesting 
an orbital period of around 25 d. 
Here we present timing and spectral analyses of RXTE 
space.Meffia^i3m®!QfadfelAii4.iHi+03 throughout the 2003 out- 
burst, as well as results rrom a search for the optical/IR 
counterpart. 
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2. OBSERVATIONS 

We made observations of 4U 1901+03 with the Pro- 
portional Counter Array fPCA: ljahoda et al.lll996|) and 
the High-Energy X -ray Timing Experiment (HEXTE; 
iGruber et all Il996|) instruments aboard RXTE. The 
PCA consists of 5 Proportional Counter Units (PCUs) 
each with a collecting area of ~ 1400 cm 2 and a 1° field 
of view, that are sensitive to X-ray photons with ener- 
gies in the range 2.5-90 keV. Photon arrival times are 
measured to « 1 /is, while spectra are accumulated in 
up to 256 energy channels. The HEXTE comprises two 
clusters, each with 4 scintillation detectors sensitive to 
photons in the range 15-250 keV, collimated to view a 
common 1° field. The detectors in the two clusters pro- 
vide a total collecting area of 1600 cm 2 . 

Short (w 3 ks) observations were scheduled every 3-4 
days between 2003 February 10 and July 16 (MJD 52,680 
and 52,837) in order to adequately sample the 25 d can- 
didate orbital period. In addition, several longer obser- 
vations were scheduled near the peak of the outburst as 
part of a separate proposal to search for cyclotron res- 
onance features (PI: Heindl). Data were analysed using 
lheasoft release 5.3 (2003 November 17). We extracted 
PCA and HEXTE spectra from intervals within each 
observation during which the center of the field-of-view 
was within 0.02° of the position of 4U 1901+03, and for 
which the limb of the Earth was more than 10° from the 
source direction. Spectra were extracted from standard 
observing mode data for each instrument ( "Standard-2" , 
with 129 channels between 2-60 keV for the PCA, and 
"Archive" , with 64 channels between 15-250 keV for the 
HEXTE). PCA spectra were accumulated separately for 
each PCU, and instrument response matrices were gen- 
erated for each PCU and each observation using PCARSP 
v.10.1. We estimated background count spectra for the 
PCA using "bright" source models (suitable for when the 
count rate exceeds « 40 counts s _1 PCU -1 ) developed for 
gain epoch 5 (2000 May 13 onwards) with PCABACKEST. 
We measured the mean flux for each observation by fit- 
ting spectra from individual PCUs separately (using the 
model described in A3. If) between 2.5-25 keV 2 . We cor- 
rected the integrated 2.5-25 keV flux (except for PCU 2) 
by dividing by the mean ratio of the fluxes for each PCU 
relative to PCU 2, and adopted the residual standard de- 
viation on the rescaled fluxes as the la uncertainty. For 
our timing analysis we used full-range PCA lightcurves 
with 4 ms time resolution, with time bins corrected to 
the solar system barycenter. 



3. RESULTS 

3.1. Flux evolution and X-ray spectrum 

The X-ray flux peaked at almost 8 x 10 -9 ergs cm" 2 s _1 
(2.5-25 keV; approximately 240 mCrab) around 2003 
February 20 (MJD 52,690) and decreased linearly down 
to < 10" 10 ergscm^s- 1 by 2003 July 15 (Fig. [TJ). Al- 
though bright outbursts in transient pulsars are some- 
times followed by extended periods of lower-lev el ac- 
tivity (e.g. KS 1947+300: [Galloway et all l2004|) . the 
ASM did not detect 4U 1901+03 at a significant 

2 The Crab flux in this energy range is 3.3 X 10 — 8 ergs cm -2 s — 1 
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Fig. 1— X-ray intensity of 4U 1901+03 throughout the 2003 
outburst, as measured by the RXTE/PCA. Error bars indicate the 
la uncertainties. 

level between the end of the 2003 outburst through 2005 
June. The modest PCA energy resolution at low en- 
ergies meant that it was not possible to constrain the 
column density nn for neutral absorption in the spec- 
tral fits. Thus, for all our spectral fits we froze nn 
at 1.2 x 10 22 cm -2 (at the lower end of the expected 
range inferred from the Ay estimates; see i)3.4|l . Our 
broadband (absorbed) flux measurements were not sen- 
sitive to the assumed value of n#. Fits with commonly- 
used pulsar X-ray spectral models including power law, 
cutoff power law and a combination of blackbody and 
power law gave x 2 values indicating statistically unac- 
ceptable fits. We found the best agreement with the data 
for a model consis ting of a Com ptonisation component 
(comptt in XSPEC: lTitarchulll994[l and a Gaussian com- 
ponent centered around 6.4 keV to represent fluorescent 
Fe line emission, both attenuated by neutral absorption 
with column density at the survey value. We assumed 
a systematic error of 1% in order to achieve a reduced 
X 2 (averaged over the fits to spectra from each PCU) 
of » 1. A systematic error of this magnitude is typi- 
cally required for PCA spectral fits to bright sources, for 
example the Crab pulsar (R. Remillard, pers. comm.). 
Near the end of the outburst (between MJD 52,780 and 
52,820) the reduced-x 2 was somewhat larger, between 2 
and 3.5. The main factor contributing to the poor x 2 
value was a deficit of photons (compared to the model) 
between 8 and 10 keV; this deficit was observed irre- 
spective of the choice of continuum components tested. 
In low signal-to-noise spectra (e.g. from short observa- 
tions), the residuals could be removed by including a 
blackbody component with fcTbb ~ 1 keV, but for longer 
observations the residuals were more complex. 

Although the residuals indicate that the adopted spec- 
tral model does not completely describe the source spec- 
trum, the derived parameters present a qualitative de- 
scription that is adequate for flux measurements as well 
as a general characterisation of the spectral shape and its 
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Fig. 2. — Representative photon spectra of 4U 1901+03 from the 
peak (February 16, or MJD 52,686) and late in the decay (June 28, 
MJD 52,818) of the 2003 outburst. The effect of the instrumental 
responses have been removed ("unfolded"). Data from PCU #2 
are shown, along with data from HEXTE clusters 1 and 2. For 
the 2003 June 28 observation, we rebinned the HEXTE spectra by 
a factor of 4. Error bars indicate the la uncertainties. The top 
and bottom panels show the residuals to each model fit (units of 
normalized counts s -1 keV - 1 ) for the 2003 February 16 and 2003 
June 28 spectra, respectively. 

variation. The fitted optical depth r decreased slightly- 
over the course of the outburst, from 6 at the peak to 
4 near the end of the outburst. The temperature of 
the scattering electrons kT e increased over the same pe- 
riod from 4 to 7 keV. A narrow emission feature around 
6.4 keV was present throughout, with equivalent width 
between 60 and 150 eV. 

We also made combined fits to PCA and HEXTE spec- 
tra for selected observations near the peak and end of the 
outburst. The spectral parameters for these broadband 
fits were similar to the fits to the PCA data only. The 
spectrum was rather soft for a HMXB pulsar, with lit- 
tle emission detected above 80 keV. As indicated by the 
evolution of the Comptonisation model parameters, the 
spectrum hardened considerably over the course of the 
outburst (Fig. EJ . Although the best-fit residuals still in- 
dicated systematic deviations at energies < 10 keV from 
the best-fit model spectrum, we found no evidence for 
cyclotron resonance features in the spectrum. Using the 
broadband fits, we estimated the bolometric correction 
as the ratio of fluxes integrated over an idealised response 
matrix spanning 0.1-200 keV, and the flux in the range 
2.5-25 keV, as 1.12. 

3.2. Pulse timing 

We estimated the pulse frequency for each observation 
by first folding the 4-ms lightcurve on a trial period, to 
obtain an observation-averaged pulse profile with (typi- 
cally) 32 phase bins. We then folded individually 256-s 



segments on the same period, and cross-correlated the re- 
sulting pulse profiles with the observation-averaged pro- 
file to obtain the phase delay for each segment. We then 
adjusted the period and repeated the procedure until the 
phase delay exhibited no net trend with time throughout 
the observation. The error was estimated from the uncer- 
tainty on the first-order term of a linear fit to the phase 
delays. The resulting frequency history shows approxi- 
mately sinusoidal variations indicative of Doppler shifts 
from binary orbital motion, superimposed on a signif- 
icant (non-linear) spin-up trend over the course of the 
outburst (Fig. EJ). 

We fit the frequency measurements with a linear model 
comprising the spin-up due to accretion torques in ad- 
dition to the apparent changes due to orbital Doppler 
shifts: 



f{t) — fspin (t) 



2irfoax sini 



forb 

x (cos I + g sin 21 + h cos 21) 



(1) 



where f sp i n (t) is the time-dependent neutron-star spin 
frequency, /o is a constant approximating f sp in(t), 
ax sini is the projected orbital semimajor axis in units 
of light travel time, and P or b is the orbital period. The 
coefficients g (— esinw) and h (= ecosw) are functions 
of the eccentricity e and the longitude of periastron u>. 
Finally, I = 2-n{t — T w / 2 )/Porb + tt/2 is the mean longi- 
tude, with Ttt/2 the epoch at which the mean longitude 
I = 7r/2. For a circular orbit, 7^/2 is the epoch of superior 
conjunction (when the neutron star is behind the com- 
panion). The right-most term in Eqn. 2 represents the 
orbital Doppler shifts to first order in e; given the mag- 
nitudes of the uncertainties of our measurements, this 
should be an adequate approximation as long as e < 0.2. 

We described the intrinsic spin frequency evolution 
with both / and / terms: 



fs P in(t) = fo + f{t-t ) + f(t-t f 



(2) 
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Fig. 3.— Pulse frequency of 4U 1901+03 throughout the 2003 
outburst. The top panel shows the measured frequency for each 
observation (open squares) with error bars indicating the la un- 
certainties. The best-fit orbital model (Table is plotted as a 
solid line; the model for the intrinsic spin evolution of the pulsar 
(excluding the orbital model) is plotted as a dashed line. The lower 
panel shows the residuals to the orbital model. 
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Fig. 4. — Pulse arrival time delays due to the 22.6 d binary orbit, 
with respect to a model for the intrinsic spin evolution of the pulsar. 
Squares indicate the measured time delays, the solid curve is the 
best-fit orbital model, and the diamonds show the residuals from 
the total model (on a 15 X expanded vertical scale). For a circular 
orbit, the mean longitude is the orbital phase angle measured from 
the ascending node; the pulsar is behind the companion at 90°. 
For such a low-eccentricity orbit, this remains approximately true. 

where fo and to are the frequency and time, respectively, 
of the first frequency measurement and the /, / are con- 
stant over the outburst duration. 

We fit the frequency model to the measurements us- 
ing a nonlinear gradient-expansion algorithm (curvefit 
in IDL). We achieved an acceptable fit (x 2 = 77.22 for 
53 degrees of freedom) with orbital parameters P or b = 
22.58 d, a x sini = 106.9 lt-sec and e = 0.035. We used 
this preliminary orbital model to perform a fit of the 
accumulated phase delay measurements. We first inte- 
grated equations n] an d [21 to obtain an expression for the 
pulse phase evolution with time: 

4>(t)=4>o + fo(t - t ) + |(t- t ) 2 + l(t- t ) 3 



-fo a x sinz (sin I 



§ cos 21 



i2l) (3) 



where 4>o is an arbitrary reference phase, corresponding 
in this case to the peak of the fundamental of the first 
pulse observed in the first observation. We then fit this 
model to the measured pulse arrival times, defined as the 
peak of the fundamental Fourier component. 

The frequency of observations during the first full or- 
bital cycle of the outburst (between MJD 52,680 and 
52,700) was such that we were able to unambiguously 
track the pulse phase over the entire cycle. Even so, we 
found significant residual phase delays which varied sys- 
tematically on a timescale of a few days, with an rms 
amplitude of w 0.1 cycles. Although variations in the 
pulse profile over the entire outburst (see H3.3|) may con- 
tribute to the residuals to the phase fit, residuals were 
also present during intervals when the pulse profile shape 
was relatively consistent. Thus, there is significant in- 
trinsic timing noise present, perhaps arising from small 
changes in the instantaneous accretion rate. 



Beyond MJD 52700 the 2-3 d gaps between the RXTE 
observations introduced the possibility of pulse count 
ambiguities, although only of magnitude ±1 cycle in gen- 
eral. We repeatedly computed the phase fit after adding 
or subtracting a cycle within the data gaps, in order to 
minimise the total \ 2 until no further improvement was 
possible. Because of the timing noise, the resulting % 2 
calculated using the errors on the individual phase mea- 
surements was much larger than the number of degrees 
of freedom. In order to estimate the confidence limits for 
the orbital parameters, we re-scaled the pulse arrival time 
errors so that the resulting x 2 was 1 per degree of free- 
dom. We then varied each parameter in turn, fitting with 
all other parameters free to vary, to determine the param- 
eter range for which the rescaled \ 2 < Xmin(l + 1/^), 
where n is the number of degrees of freedom (1473 for 
the full set of arrival time measurements). The resulting 
orbital parameters and uncertainties are listed in Table 
the predicted frequency, intrinsic spin frequency for 
the neutron star and the residuals from the model are 
shown in Fig. [3] The pulse arrival time delays with re- 
spect to the intrinsic spin evolution model are shown in 
Fig. H 

Our best-fit parameters describing the intrinsic spin 
evolution indicate that the spin-up rate was initially 
around 3 x 10~ n Hzs -1 , similar to the maximum mea- 
sured for other transient pulsars (e.g. iBildsten et alJ 
Il997|) . According to the intrinsic spin model, the spin- 
up decreased throughout the outburst, falling to zero just 
before the end of the outburst, around MJD 52,824. The 
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Fig. 5. — Representative pulse profiles throughout the 2003 out- 
burst of 4U 1901+03. Phase 0.0 corresponds to the maximum of 
the fundamental. The profiles have been rescaled by a factor of 3 
to exaggerate the pulse amplitude. 
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2.5-25 keV flux by this time had dropped to around 5.5 x 
10 -10 ergscm _2 s _1 . 

3.3. Pulse profile variability 

The pulse profile was consistently non-sinusoidal, and 
exhibited several intervals of stability punctuated by rela- 
tively rapid change. The fractional pulse amplitude also 
varied over the outburst, and was generally between 4 
and 22% (rms). From the beginning of the outburst un- 
til 2003 March 3 the profile was consistently similar to 
the example shown from February 19 (MJD 52,689; Fig. 
|SJ). During that observation the pulse amplitude was 15% 
rms. Between March 3 and 13 the profile switched to the 
characteristic double-peaked shape of March 19; also on 
March 13 the rms amplitude dropped to 7%. The profile 
continued to evolve throughout March and April, with 
gradually increasing pulse fraction and small drifts of the 
harmonic components until another abrupt shift between 
2003 June 14-18. The pulse amplitude had risen to 13% 
rms on June 14 before falling abruptly to 9% rms on June 
18, and then recovering to an overall maximum for the 
outburst of 22% rms on June 28. After this final peak, 
the profile remained similar in shape but with steadily 
decreasing amplitude towards the end of the outburst. 
Pulsations became undetectable (< 1% rms) after 2003 
July 13 (MJD 52,833), by which time the 2.5-25 keV flux 
had dropped to below 3 x 10~ n ergscm _2 s _1 . 

3.4. A search for the optical/IR counterpart 

Following th e improved X-ray posit ion obtained from 
the PCA scan (Gallowa y et alJ[2003aj) . we examined the 
field of 4U 1901+03 in Digital Sky Survey and 2MASS 
images to identify the optical counterpart. The mass 
function implies a minimum companion mass (assuming 
a 1.4M© neutron star) of 4.5M Q . The most probable 
companion mass (for an inclination of 60°) is 6.0 M . 
The donor stars in long-period binary pulsars are typi- 
cally either Be stars or (sometimes Roche-lobe filling) OB 
supergiants. The column density interpolated from Hi 
survey observation s towards the source is between (1.1- 
1.24) x 10 22 cm" 2 ■: Dukev fc Lockmanlll990HStark et"aTl 
Il992|) . which translates to an Ay = 6.1-6.9 (assum- 
ing a standard dust to gas ratio; iPredehl fc SchmitH 
1995). Alternatively, the reddening estimated from dus t 
IR emission is higher at A v = 10.2 iSchlegel et alJ[l 998). 
The USNO A2.0 astrometric catalog ijMonet et al.lll998j) 
contains just two stars within the 1' error circle of 
4U 1901+03 having colors consistent with early-type 
stars suffering extinction with Ay > 6 (i.e. B — R > 3; 
Fig. [BJ upper panel). We obtained low-resolution spectra 
of these two candidates using the Low Dispersion Survey 
Spectrograph (LDSS-2) on the 6.5 m Clay (Magellan II) 
telescope at Las Campanas, Chile. We accumulated two 
600 s spectra of each candidate on 2003 August 10, using 
the medium red grism with a 1" long slit, covering the 
range 4500-9000 A. However, the overall spectral fea- 
tures suggest that these candidates are instead low-mass 
K0-5 stars (Fig. H3 lower panel), which are ruled out as 
the counterpart on the basis of the X-ray mass function. 
Clearly, it is not possible to distinguish between early- 
type stars with high extinction and nearby late- type stars 
from optical photometry alone. 

Thus, we also examined the J, H and K magnitudes 
of candidates within 1' of the X-ray position from the 



2MASS point source catalogue 3 . The expected range of 
colors for a B-star for the estimated exti nction ran ge are 
J-H = 0.5-1.1 and H - K = 0.3-0.6 lICo^EOQO). The 
colors expected for a supergiant counterpart fall within 
similar ranges, slightly larger in J — H and smaller in 
H — K. We found approximately ten stars in the 1' 
RXTE error circle with colors within these ranges, in- 
cluding the optically identified candidates A and B. Star 
B was the brightest in the IR bands, with J = 11.5; this 
is the only star consistent with the expected brightness 
of a supergiant counterpart at < 10 kpc, and can be 
ruled out as the counterpart on the basis of our spec- 
troscopic observations. Two other stars had J mag- 
nitudes similar to that of star A, at J fa 13; we ex- 
pect these candidates are also nearby low-mass stars, 
like star A. The remaining candidates all had J > 14, 
much fainter than the limit of J ~ 12 expected for 
a supergiant companion at d < 10 kpc. Several of 
these stars were not detected in the DSS image, sug- 
gesting that B > 20 perhaps consistent with the up- 
per end of the estimated Ay range towards the source. 
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Fig. 6.— Upper panel Field of 4U 1901+03 from the DSS POSS- 
II data, showing the 1' RXTE error circle and the two candidate 
counterparts A and B identified by their large B — R color indices. 
The B-magnitudes measured from the USNO-A2.0 catalog, are 
18.5 and 19.3, respectively; star B is also the brightest star within 
the 1' error circle in IR at J = 11.5. Lower panel Magellan spectra 
of the candidates A a nd B, as well as a comparis on spectrum from a 
K2V star, HD 22049 IDanks fc Dennefeldll99l) . The spectra of A 
and B were dcrcddencd with E(B — V) = 1.1 and 1.5 respectively 
IFitzpa trick 1999|). The spectra of stars A and B are consistent 
with KO-5 stars and thus can both be ruled out as the counterpart 
of 4U 1901+03. 

3 \protect ittp://www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.h 
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4. DISCUSSION 

4U 1901+03 is one of a small, but growing class of low- 
eccentricity high-mass X-ray binaries. Although we were 
unable to identify the optical counterpart, the limit of 
J > 13 for stars in the 1' RXTE error circle rules out a 
supergiant companion, unless the distance to the source 
is > 10 kpc. 

Furt hermore, the source is located on the Corbet dia- 
gram l)Corbetll986D with confirmed Be transients such as 
4U 0115+63 fe.g. lBildsten et alJll997|) . Like that source, 
4U 1901+03 exhibits infrequent outbursts that can span 
multiple orbital periods, a behavior quite unlike the typi- 
cally persistent activity of partially Roche-lobe filling su- 
pergiants. If the mass donor is not filling it's Roche lobe, 
the efficacy of tidal forces in circularizing such a wide 
orbit is negligible, which makes the present low eccen- 
tricity of 0.036 difficult to understand (given the initially 
eccentric orbit expected to arise as a result of the natal 
supernova kick). As suggested bv lPfahl et alJ l)2002|) . the 
formation events for low-eccentricity O-B transients like 
4U 1901+03 may be dynamically distinct from the more 
common high-eccentricity binaries due to a much lower 
initial kick to the neutron star. 

The lack of eclipses indicates the inclination is < 
85°. The 95% upper limit on the companion mass for 
an a priori isotropic distribution of inclination angles 
is 88 Mq. The estimated peak bolometric luminos- 
ity (for a distance of 10 kpc) was 1.1 x 10 38 ergss -1 , 
and the integrated luminosity over the course of the 
outburst was 7.4 x 10 44 (<i/10 kpc) 2 ergs (for a neutron 
star with R = 10 km and M = 1.4 Mq). The max- 
imum intensity observed by Uhuru during the 1970-1 
outburst was 87+ 11 counts" 1 at epoch 1971.0, cor- 
respondin g to 1.5 x 10 -9 ergs cm" 2 s" 1 in the range 
2-6 keV (|Forman et al.lll976h . Our broadband RXTE 
spectra suggest that 20-30% of the source flux is emit- 
ted in this energy range, so that the estimated maxi- 
mum bolometric luminosity for the Uhuru observations 
was w 8 x 10 37 ergss -1 (for d — 10 kpc), consis- 
tent with the peak RXTE value to within the error. 
Vela 5 B measured a peak of 8 cou nt s" 1 around epoch 
1970.9 (|Priedhorskv & Terrel]lll984jK The maximum ob- 
served flux was thus 3.6 x 10 -9 ergs cm" 2 s" 1 in the 
range 3-12 keV, in which range 60% of the flux ob- 
served by RXTE is emitted; thus, the estimated max- 
imum bolometric luminosity measured by Vela 5B was 
7 x 10 37 ergs cm" 2 s" 1 , roughly consistent with both 
the Uhuru measurements of the 1970-1 outburst and 
the the peak measured by RXTE for the 2003 out- 
burst. We also note that the estimated duration of 
the 1970-1 outburst (excluding the peak around epoch 
1970.6) was at least 120 d, similar to the 150 d dura- 
tion of the 2003 outburst. No emission prior to 2003 
February was detected by the RXTE/ ASM, although 
there is the possibility that one or more intervening out- 
bursts occurred sometime between 1971 and the launch 
of RXTE in 1995. A similarly bright outburst before 
1980 would probably have been detec ted by Vela 5B or 
the ASM onboard Ariel 5 l|Holtlll976|) . An outburst be- 
tween 1987 and 1996 would likely have been detected 
with the Ginqa ASM (operational until 1991 November; 
iTsunemi et al.lfl989^ or the BAT SE experiment onb oard 
CGRO (1991 April to 2000 ,Tune: IZhang et alll995ft . As- 



suming no intermediate outbursts occurred, we derive a 
time-averaged accretion rate for a 32.2 yr recurrence time 
of 8.1 x 10" n (d/10kpc) 2 M Q yr -1 . We note that the 
32.2 yr interval between outbursts in 4U 1901+03 may 
be the longest presently known for any X-ray transient. 

Whilst no cyclotron absorption features were detected 
in the X-ray spectrum, the observed range of source flux 
over which pulsations were detected allows a rough esti- 
mate of the dipole magnetic field strength of the neutron 
star. Assuming that a disk is present, it must be trun- 
cated above the surface of the neutron star in order to 
allow the magnetic field to channel the accreting mate- 
rial and produce observable pulsations. The truncation 
radius T m is inversely propo rtional to the mass accretion 
rate (e.g. lFrank et alJll9 92'l. so that the requirement for 
Tm > i?» (where i?» ks 10 km is the neutron star radius) 
even at the peak of the outburst implies a lower limit 
on the magnetic field strength (although for 4U 1901+03 
this limit is several orders of magnitude below the canon- 
ical field strength for long-period pulsars of ~ 10 12 G). 
For accretion to be dynamically feasible also requires that 
the inner disk radius is within the corotation radius r co , 
i.e. the radius at which the Keplerian orbital frequency 
equals the neutron star spin frequency. For the lowest 
flux at which pulsations were detected, this implies that 
B < 0.5 x 10 12 (d/10kpc) 2 G, giving a fundamental cy- 
clotron frequency of v cyc < 4(d/10kpc) 2 keV. This is 
consistent with the absence of detectable cyclotron fea- 
tures in the broadband spectrum, although it is possi- 
ble that the residuals frequently present at < 10 keV 
may arise from higher harmonics of a low-energy cy- 
clotron absorption line. We can also estimate the mag- 
netic moment by assuming that the long-term mass ac- 
cretion has left the pulsar close to s pin equilibrium, i.e. 
tm ~ r co (e.g. iBildsten et al.l ll997'). In that case, we 
find B w 0.3 x 10 12 (d/10kpc)- 6 / 7 G, consistent with 
the above estimate. We note that the inferred limit on 
v cyc and the low values of kT e are qualitatively consis- 
tent with the observed correlation in other pulsars be- 
tween the cyclotron energy and the spectral cutoff (e.g. 
iCoburn et al]l200^) . 

Identification of the counterpart to 4U 1901+03 is es- 
sential to confirm the nature of the mass donor suggested 
by these observations, which may prove difficult unless 
the source becomes active again. If the outbursts occur 
regularly every « 30 yr, significantly more advanced X- 
ray instruments may be available for the next outburst, 
perhaps sufficient to measure the position more precisely 
and identify the optical counterpart, as well as resolve 
the residuals below 10 keV and measure the properties 
of the low-energy cyclotron lines, thus allowing direct 
measurement of the field strength. 

This research has made use of data obtained through 
the High Energy Astrophysics Science Archive Research 
Center Online Service, provided by the NASA/Goddard 
Space Flight Center. The Second Palomar Observa- 
tory Sky Survey (POSS-II) was made by the California 
Institute of Technology with funds from the National 
Science Foundation, the National Geographic Society, 
the Sloan Foundation, the Samuel Oschin Foundation, 
and the Eastman Kodak Corporation. This publication 
makes use of data products from the Two Micron All 
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TABLE 1 

Orbital parameters from the pulse arrival time fit for. 
4U 1901+03 



Parameter 


Value a 


^orb (d) 


22.5827 ± 0.0002 


ax sini (lt-s) 


106.989 ± 0.015 


T n/2 (MJD) 


52682.2440 ± 0.0010 


e 


0.0363 ± 0.0003 


^peri ( ) 


268.812 ± 0.003 


MM) (M ) 


2.58348 ± 0.00005 


to (MJD) 


52680.33163 


fo (mHz) 


361.970908 ± 0.000009 


/ (10- 11 Hz s- 1 ) 


2.9759 ± 0.0004 


/ (10- 18 Hz s-2) 


-1.1997 ±0.0003 



a The uncertainties are at the la confidence level, i.e. the range of 
each parameter for which the rescaled \ 2 < Xmin + ^Ln/^P ^)- 



Sky Survey, which is a joint project of the University of 
Massachusetts and the Infrared Processing and Analysis 
Center /California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the 



National Science Foundation. This work was supported 
in part by the NASA Long Term Space Astrophysics pro- 
gram under grant NAG 5-9184. 
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